Statistical modelling of electrochemical 

■^- . deposition of nanostructured hybrid films 

§ : with ZnO-Eosin Y as a case example 

(N 

^: 

< ; F. D. A. Aarao Reis^'* J. P. Badiali^ Th. Pauporte' 

g;j: D. Lincot^ 



(N 
> 



' Instituto de Fisica, Universidade Federal Fluminense, Avenida Litordnea s/n, 

24210-340 Niteroi RJ, Brazil 

^Laboratoire d'Electrochimie et de Chimie Analytique (UMR 1515 - CNRS - 
ENSCP - Paris 6), Universite P. et M. Curie, 4 Place Jussieu, 15005 Paris, 



a 

CD 

r~| ■ France 

• • ^ Lahoratoire d'Electrochimie et de Chimie Analytique (UMR 1515 - CNRS - 

Q , ENSCP - Paris 6), Ecole Nationale Superieure de Chimie de Paris, 11 rue Pierre 
' c/3 ' et Marie Curie, 15231, Paris, France 



Abstract 



We study models of electrodeposition of hybrid organic-inorganic films with a special 

focus on the growth of ZnO with eosin Y. First we propose a rate equation model 

Cp . which assumes that the organic additives form branches with an exposed part above 

the ZnO deposit, growing with larger rate than the pure film and producing ZnO 

\^ ' at the exposed length. This accounts for the generation of OH~ ions from reduction 

^^ ■ of dissolved oxygen near the branches and reactions with Zn ~^ ions to form ZnO 

y^ , molecules. The film grows with the same rate of the branches, which qualitatively 

• 1—1 I explains their catalytic effect, and we discuss the role of the additive concentration. 

J>^' Subsequently, we propose a statistical model which represents the diffusion of the 

r^ . hydroxide precursor and of eosin in solution and adopt simple probabilistic rules for 

'"^' the reactions, similarly to diffusion-limited aggregation models. The catalytic effect 

is represented by the preferencial production of OH~ ions near eosin. The model 

S^ . is simulated with relative concentrations in solution near the experimental values. 

H I An improvement of the growth rate is possible only with a rather large apparent 

diffusion coefficient of eosin in solution compared to that of hydroxide precursors. 

When neighboring eosin clusters competitively grow, the increase in the growth rate 

and a high eosin loading are observed in the simulated deposits. Those features are 

in qualitative agreement with experimental results. 

Key words: Zinc oxide, Eosin Y, cathodic electrodeposition, growth model, 
diffusion-limited aggregation 
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1 Introduction 



The growth of thin films from solutions has emerged as an efficient, low tem- 
perature, versatile preparation route which can be used at large scale for the 
production of high quality materials [1]. One of the interests of these methods 
is the possibility of adding foreign soluble compounds to the deposition bath 
and modulating film properties by playing with the interactions which arise 
between the growing film and these compounds. Among the additives, organic 
ones are of utmost interest since many of them have been shown to act as 
templating agents and/or crystal growth directing agents, and in many cases 
as functionalizing agents for the deposits. 

These effects are well-documented in the case of zinc oxide. This inorganic 
compound has attracted much attention due to a broad range of potential high 
technology applications such as surface acoustic wave filter [2], light emitting 
diodes [3], lasers [4], varistors [5], gas sensors [6] and solar cells [7]. ZnO can 
be prepared as a thin film by chemical deposition methods [8,9] or by elec- 
trodeposition [10,11,12,13,14,15]. The dramatic effects of organic additives on 
film structures and morphologies have been reported with both deposition 
routes. For instance, in the presence of citrate, chemically deposited ZnO are 
formed of stacked nanoplates assembled in a biomimetic manner [9]. However 
the most impressive morphological changes have been observed in electrode- 
posited zinc oxide. In the presence of sodium dodecylsulfate (SDS), the for- 
mation of lamellar nanostructures has been shown [16,17]. The SDS can be 
removed from the film and highly porous thin films are released [17]. Different 
dyes, such as eosin Y (EY) [18,19,20,21], fluorescein (FL) [22], Tetrasulfonated 
phthalocyanines (TSPc) [23,24,25,26] and riboflavin (RI) [27], have also been 
shown to incorporate in the fllm and to give rise to nanostructures. These 
organic dyes contain negatively charged functions (carboxylate (EY, FL), sul- 
fonate (TSPc) or phosphonate (RI) groups) which allow their direct binding 
to the oxide crystal surface during the synthesis. 

While the effects of additives in fllm shaping are widely illustrated in exper- 
imental works, particularly that of the organic ones, the exact role of these 
compounds in the growth process remains to be clarifled in most cases. The 
aim of the present work is to flll this gap by proposing kinetic and statisti- 
cal models for zinc oxide electrodeposition in the presence of EY in solution, 
which focus on a small number of basic features of those processes. The start- 
ing point for this modeling is a series of experimental results, some of them 
also observed in deposition with different organic additives [17]. 
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The first step of electro diemical zinc oxide syntliesis is tlie catliodic reduction 
of a hydroxide precursor such as molecular oxygen [11,13], hydrogen peroxide 
[14,15] or nitrate ions [10,18]: 

]^02 + H20 + 2e- ^20H-, (1) 

H2O2 + 2e- -^ 20H~ (2) 



or 



N0^~ + H20 + 2e~ ^ NO2' + 20H'. (3) 



At a temperature slightly above room temperature [28], well crystallized zinc 
oxide with the wurtzite hexagonal structure is precipitated at the electrode 
surface: 

Zn^+ + 20H- -^ ZnO + H2O. (4) 



In many cases, it has also been shown that the organic component presents 
catalytic properties for the deposition process. This catalytic effect is con- 
nected to the improvement of O2 reduction (Eq. 1) with EY [20], TSPc [26], 
SDS [17], whereas EY has also been shown to catalyze the H2O2 reduction 
(Eq. 2) [19]. Thus, one of the assumptions of our models is a distinguished 
rate of formation of ZnO near the species that represent aggregated eosin 
molecules. 

The zinc oxide films prepared in the presence of EY show large round shaped 
single crystals of zinc oxide filled with self-assembled dye aggregates, as shown 
in Fig. 1 [21]. It is observed that the dye can be almost completely removed 
by a soft chemical treatment [29,30], revealing a network of mesopores which 
is connected to the surface and can be filled with a solution. This network, 
formed by aggregated dye molecules, acts as a template for the growth of ZnO. 
In our models, these experimental facts justify the assumption of formation 
of connected branches of the species that represent the additive molecules in 
the films. 

In Sec. II of this paper, we will introduce our first model, which is based 
on rate equations, and explains some features of the cooperative growth of a 
ZnO film and additive aggregates. This model is particularly useful for un- 
derstanding the role of additive concentrations and diffusion coefficients in a 
qualitative way. Subsequently, in Sec. Ill we will present a statistical model 
which represents the main microscopic features of those processes by adopting 
probabilistic rules of diffusion and aggregation of some chemical species in 



solution and in the deposit. The hypothesis of a growth mechanism controlled 
by diffusion-limited aggregation of the eosin molecules was anticipated in Ref. 
[31]. Simulation results of that statistical model, presented in Sec. IV, provide 
estimates of growth rates, structures of ZnO and EY deposits and shows a 
role of diffusion coefficients which are in qualitative agreement with experi- 
mental findings. Finally, in Sec. V we summarize our results and present our 
conclusions. 



2 Rate equation model 



This model aims at explaining basic kinetic features of the growth of ZnO 
with organic additives with drastic simplifications of the film structure. 

First we assume that pure ZnO growth takes place with rate tq (in nanometers 
per second), due to reactions of type (4) near the film surface. If the average 
volume occupied by one molecule is V , then tq/V is the growth rate per unit 
area in that case. Even in the presence of the additives, the roughness of the 
ZnO film surface will be neglected. 

When the additive is present, we assume that it forms straight structures 
which grow from the electrode, as shown in Fig. 2. This is consistent with the 
experimental finding of porous structures when eosin Y was removed from the 
hybrid films [29,30], but with an initial oversimplification of that structure. 
The total number of linear branches of the additive is Nq, which depends on 
the nucleation of the additive at the electrode surface in the beginning of the 
growth process. Each branch is assumed to grow with rate ta if it is not covered 
by the ZnO deposit (r^ > tq), so that the height of each branch increases as 
H = rxt. ta is expected to increase with the concentration p^ and the diffusion 
coefficient Da of the additive in solution, but the exact dependence on those 
quantities will not be important at this point. 

In order to represent the catalytic effect of the additive, we assume that an 
excess of OH is produced near the branches, which leads to the production 
of ZnO at those regions. The rate of production of ZnO is expected to be 
proportional to the length of those branches which is exposed above the film 
surface, {H — h) (see Fig. 2), where h is the ZnO film height. Consequently, 
the total rate of production of ZnO molecules near the additive branches 
is Nok{H — h), where /c is a reaction constant. A constant reaction rate in 
this model accounts for the assumption that the oxygen reduction reaction 
is not under diffusion control. Also, it is assumed that the ZnO molecules 
immediately precipitate and aggregate to the flat film surface, leaving the 
branches uncovered. 



Under these conditions, if the area of the film is A, then the number of ZnO 
molecules deposited per unit time is 

^ = yA + Nok{H-h). (5) 



Here, the first contribution in the right hand side comes from the reactions 
ocurring near the film surface and the second one comes from the reactions 
near the branches. The growth rate of the ZnO deposit is ^ = ^^, which 
gives 

dh , NokV 

— =ro- ah + TAat , a = — - — . (6) 



The solution of Eq. (6) is 

/. = r^t + i!:A^I^(e---l). (7) 



It gives a constant growth rate r^ for the deposit at long times, with 1/a being 
the characteristic time of decay. It means that the film growth rate attains the 
same growth rate of the branches in the steady state. The rate of production 
of ZnO near the film surface is ro < r^, but the exposed height of the branches 
saturates aX H — h = ''^~''° , and produces an excess of ZnO molecules in their 
neighborhood which is enough for branches and film to grow with the same 
rate. 

The dependence of Vj^ on A'"o is not explicit in the model, but it exists. As 
the number of branches per unit area increases, the aggregation of the dif- 
fusing additive to them is facilitated and, consequently, their growth rate ta 
increases. From the dependence of the saturation value of H — h on a, we ex- 
pect that it decreases as A'"o increases: the steady state regime is more rapidly 
attained (1/a decreases), with a smaller exposed height of the branches. On 
the other hand, the dependence of if — /i on r^ suggests the opposite behavior, 
i. e. that it increases with Nq. The interplay of these effects will determine the 
overall dependence oi H — h on Nq. Anyway, this analysis shows the nontrivial 
consequences of the nucleation processes (which determine Nq) on the long 
time behavior of the system. 

Now let us consider the effect of the additive concentration in solution, p^, 
assuming that its diffusion coefficient is constant (constant temperature con- 
ditions). Since tq only represents the reactions occurring near the film surface, 
it is expected to be constant as pA increases. On the other hand, the growth 
rate of the branches, ta, is some monotonically increasing function of pA- For 
low concentrations, ta is small, consequently the branches will be covered 



by ZnO. Quantitatively, this is the case where ta [pa) < f^o, which imphes 
Pa < Pc, where pc is a critical concentration of additive such that r^ (pc) = ''"o- 
For higher concentrations, we have r^ > Tq and the above solution of the rate 
equation model is valid. 

The dependence of the growth rate on the additive concentration, with the 
catalytic effect only for r^ > ro, is illustrated in Fig. 3. There, we assume 
that ta linearly increases with pA when the concentration is not too large 
(first order reaction). This is the regime where our model is expected to work. 
For large concentrations. Fig. 3 shows a saturation of growth rate, which is a 
consequence of the limitations to oxigen diffusion towards the surface. Indeed, 
in this regime, the catalytic effect of eosin to reduce oxygen does not result in 
an increased production of hydroxide ions anymore. The hydroxide precursor 
current towards the electrode becomes constant, which also occurs with the 
growth rate of ZnO, which is directly linked to the availability of hydroxy de 
ions. 

The behavior observed in the intermediate EY concentration range is in qual- 
itative agreement with the findings of Ref. [19] for ZnO/EY electrodeposi- 
tion with hydrogen peroxide as hydroxide precursor: no significant effect in 
the growth rate is found for EY concentrations below 2pM, but an increase 
is observed for EY concentrations above that value. On the other hand, no 
well-defined plateau in the growth rate is observed in electrodeposition with 
molecular oxygen [32], but an increase in the growth rate even for small addi- 
tive concentrations. This is probably due to nucleation features that are not 
properly represented in the present simplified model. 

The main parameters of the rate equation model are the growth rates ta and 
ro, which depend on several parameters, such as concentrations of hydroxide 
precursors and additives in solution, zinc ions concentration and the respective 
diffusion coefficients. However, the model leads to two important conclusions. 
First, in order that the catalytic effect is observed, it is necessary that ta > tq, 
consequently the concentration and the diffusion coefficient of the additive in 
solution must be combined in a suitable way to fullfill this condition. Secondly, 
if the additive branches are able to grow above the film surface, then the film 
growth rate is the same of those branches in the stationary state (t ^ 1/a in 
Eq. 7), independently of the rate of pure ZnO films growth. These conclusions 
will guide the simulation work on the microscopic model. 



3 Microscopic model 



Now we present a more complete model for the electrodeposition of ZnO 
with organic additives, in which the physical properties of different species 



in solution and in the deposit are taken into account. This is particularly 
important to understand the role of diffusion coefficients and to predict the 
structure of the deposit, as well as to perform quantitative comparisons with 
experimental data. At the latter point, we will basically refer to ZnO/EY 
films data available, thus the presentation of the model will be inspired by 
this growth process. 

We will consider a lattice model with two species executing random walks 
in solution and two species forming a deposit grown above a fiat electrode. 
All lengths will be given in lattice units. Due to the relative complexity of 
this model, it will be solved by simulation. A two-dimensional version will be 
discussed, which we expect to capture the main features of the real process, 
at the same time of being amenable for computation. 

The species in the solution, O and E, represent the hydroxide precursor (O2 
or H2O2) and eosin Y, respectively. Their different sizes, as shown in Fig. 4a, 
qualitatively account for the different masses and radius. Their auto-diffusion 
coefficients are Dq and De, i. e. each O (E) in solution executes Dq {De) 
random movements of unit size per unit time. These rates set the time scale 
of our model. 

The deposit will be formed by two aggregated species, named ZA and EA, 
representing respectively the ZnO and the eosin Y molecules in the film, as 
shown in Fig. 4b. For simplicity, their sizes are the same as those of the species 
in solution. 

Although the zinc and hydroxide ions play an essential role in the reactions 
leading to the formation of ZnO (Eqs. 1, 2 and 4), they will not be explicitely 
represented in our model. Instead, their effects will be described by probabil- 
ities of removing O particles from the solution and producing new aggregated 
particles, thus avoiding the complications of the explicit representation of the 
above mentioned reactions. Following this reasoning, the growth process is 
represented by the simplified reactions O -^ ZA and E -^ EA, which can 
take place with given probabilities and under conditions that are related to 
the neighborhood of the particles in solution. 

The growth process begins with a flat electrode (a line of length L) and a 
solution with concentrations po and pE of O and E particles, respectively. 

The aggregation of an E particle (reaction E -^ EA) occurs when that particle 
has a nearest neighbor of a previously aggregated EA (Fig. 5a) or when it has 
a nearest neighbor site belonging to the electrode (the initial substrate). The 
new EA particle permanently sticks at the position where it is formed. The 
contact of E particles with ZA sites never generates EA particles, which is a 
key hypothesis for explaining the formation of continuous fibers of eosin Y in 
the films. 



The reaction O -^ ZA occurs upon contact of a particle O with ZA or EA. 
When the O particle has a neighbor site ZA or a neighbor site belonging to the 
electrode (initial substrate), the reaction occurs with probability p, otherwise 
the O particle remains in the solution. When the O particle has a neighboring 
EA, that reaction occurs with probability 1. These rules are illustrated in Fig. 
5b. 

Contrary to the EA particles, ZA particles are allowed to precipitate and dif- 
fuse on the film surface after their formation. In order to represent the main 
features expected for this process and save computational time, a simple diffu- 
sion mechanism will be mimicked immediately after the ZA particle is created: 
it precipitates to the topmost empty site of its current column, searches for the 
point with the highest number of nearest neighbors ZA whithin a certain ra- 
dius R (defined below), and permanently sticks at that point. These processes 
are illustrated in Fig. 5c. 

The above stochastic rules are motivated by some experimental features of the 
ZnO I'EY films growth. The rules for formation of EA particles represent the 
initial nucleation of EY at the electrode and the subsequent formation of EY 
clusters upon reduction of the molecules that reach the film surface. Since EA 
is not formed upon contact of E with ZA, there will be no isolated EA particle 
in the ZA matrix, which accounts for the fact that almost all EY can be 
removed from the hybrid films, leaving them with a porous structure [29,30]. 
As regards the formation of ZA particles, the higher probability of occurrence 
after contact of O with EA follows from the assumption that oxygen precursors 
reduction is enhanced near the eosin clusters, which facilitates the formation 
of ZnO in those regions. This accounts for the catalytic effect of eosin Y, also 
observed experimentally. 

For the diffusion of a ZA particle on the ZnO film surface after its formation, 
we consider a surface diffusion coefficient D^, which gives a radius for searching 
the final aggregation position as i? = \/TJat, where r is the current average 
time for one layer deposition. This time is updated during the simulation, 
while Da is kept fixed. The condition to choose the final aggregation point 
follows the same ideas of the Wolf- Villain model for thin film deposition [33], 
and are reasonable for formation of a crystalline structure due to the trend of 
increasing the binding energy. 

The height of the solution above the film surface is kept constant during the 
simulation of the growth process. New particles (O or E) are left at random 
positions of the top layer of the solution immediately after an aggregation 
event (rections E — > EA, O -^ ZA), in order to maintain constant concentra- 
tions. This height was 300 lattice units in our simulations. The diffusion layer 
near the film surface, in which concentrations of O and E were reduced, was of 
order of 100 lattice units or less. These conditions are reasonable to represent 



experiments in which a rotating disk electrode was used, thus avoiding the 
growth of the diffusion layer and the subsequent decrease of the growth rate. 

We recall that the rules for aggregation of E particles are the same of the 
original diffusion-limited aggregation (DLA) model of Witten and Sander [34], 
in which growth took place from a single seed. Extensions to the case of an 
initial fiat surface were also considered by several authors [35,36,37,38]. To 
be more precise, due to the collective diffusion mechanism of E particles in 
solution, the growth kinetics of E is equivalent to that of the multiparticle 
biased DLA (MBDLA) in the condition of zero bias [39,40,41,42]. Indeed, 
MBDLA was proposed as a model for electrodeposition of CoP alloys [40,43]. 
On the other hand, the aggregation rules of ZA particles, which include surface 
relaxation, suggest Edwards- Wilkinson growth in that case [44,45]. However, 
in the cooperative growth model, a nontrivial interplay of both processes is 
observed, as shown below. 



4 Results for the microscopic model 



In all simulations presented here, we considered Do = 1, and in most cases 
the lattice lenghts are L = 512 and the diffusion coefficient of aggregated ZA 
is Da = 0.005. Some data for L = 2048 were also collected in order to search 
for possible finite-size effects, mainly when the thickest deposits were grown. 
In experiments, the maximum ratio of concentrations of eosin and hydroxide 
precursor in solution is in the range 10~^ — 10~^, thus in this model we worked 
with ratios between (no eosin) and 10~^. 

The results of the rate equation model of Sec. II suggest that the growth rate 
rA of the additive (eosin) is proportional to peDe, while the growth rate of the 
pure ZnO deposit, tq, is proportional to ppoDo- In the following, the model 
parameters will be tuned in the light of these relations and the conclusions of 
the model of Sec. II. 

First the model was simulated considering De = 1/8 (thus De/Dq ~ 0.1), 
in order to account for the much larger mass of eosin when compared to the 
hydroxide precursors. Da = 0.05 was considered in all simulations, some of 
them with po = 2 x 10"^, pe = 2 x 10""^, and others with po = 5 x 10^^, 
Pe = 5 X 10~^. Several different values of p were considered, with p > 10^^. 

In Fig. 6 we show a region of a deposit obtained with p = 10~^, pz = 5x 10~^ 
and pe = 5 X 10~^, at t = 3 x 10^. Some EA particles are present at the initial 
electrode surface, but they are rapidly covered by the ZA layer. Subsequently, 
there is no possibility for uptaking other eosin molecules in the film, then the 
final growth rate is determined by ZA aggregation. From the results obtained 



with such small diffusion coefficient of E in solution, we conclude that this is 
not the suitable condition to observe the catalytic effect of the additive. 

Subsequently, we assumed that E particles had a significantly large diffusion 
coefficient, De = 1, which is equal to the coefficient of Z. In Figs. 7a, 7b and 7c, 
we show the time evolution of a region of a deposit grown with po = 2 x 10^^, 
Pe = 2 X 10~^ and p = 10~^. Although many EA particles attach to the elec- 
trode in the beginning of the process, they frequently become covered by ZA 
particles after some time, even with a small probability of 0-ZA aggregation. 
The growth process is improved only at isolated regions where clusters of EA 
particles succeed to grow, as shown in Figs. 7a- 7c. 

When compared to the case without E particles, the average height of the 
deposit increases more than twice near the clusters of Figs. 7a-c, for the same 
time of growth. However, the large distance between the surviving branches 
leads to a very large surface roughness at long lengthscales. Moreover, this 
roughness rapidly increases in time because the mounds growing with large 
rates are separated by long valleys which slowly grow. These features are not 
observed in experiments, for which we refer to Fig. 1 as a typical example. 
We also tested smaller concentrations of po, but keeping the relative concen- 
tration between E and O fixed (10^^), obtaining similar results (growth of 
EA branches may be even more difficult for the same p). We conclude that 
larger values of De must be tested in order to find the suitable conditions to 
represent the experimental results. 

Following this reasoning, we also considered the case De = 2, and obtained 
results that qualitatively agree with experimental findings. In this case, small 
concentrations of E particles are able to improve the growth process, with the 
creation of a large number of nucleation centers in the eletrode, the formation 
of a large number of branches and a competition between them as the deposit 
grows, as well as a much smaller surface roughness. The formation of an eosin 
interconnected network inside the ZnO matrix is clear, which also agrees with 
the main features of the experimental nanostructures. 

Fig. 8 illustrates the beginning of the growth of a region of a deposit with 
Po = 5 X 10~^, pe = 5 X 10"^, Da = 0.05 and p = 0.25. The average height 
for t = 3 X 10^ (Fig. 8c) is nearly 3 times larger than the height of the film 
grown without E particles during the same time, and the eosin loadings in 
the deposits (ratio between number of EA particles and ZA particles) of Figs. 
8a-8c increase from 0.64% to 0.85%. Simulation at longer times show sthat 
these quantities attain different steady state values, so that the features of 
Fig. 8 may be interpreted as typical of the initial growth process. 

One interesting point revealed in Fig. 8 is the competition between the growing 
EA clusters. The central cluster of Figs. 8a and 8b grew slower than the 
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neighboring ones and, consequently, was covered by particles ZA in Fig. 8c, 
while the larger surviving clusters created more branches. 

In Fig. 9 we show a region of a deposit grown until t = 2 x 10^ with the same 
parameters of those in Fig. 8, except that Da = 0.01. The surface roughness 
is larger in this case due to the inhomogeneous growth of EA branches and 
smaller diffusion coefficient of aggregated particles. However, other quantities 
are not much different from those obtained with D^ = 0.05, such as the growth 
rate at steady state conditions, which is nearly 3.3 times larger than the rate 
of film grown without E particles. It clearly shows the catalytic effect of the 
eosin clusters. Another important quantity is the EA loading in the film, which 
is nearly 2%. This is a high value compared to the relative concentration of 
1% in solution. 

The above results show that a realistic description of the growth of ZnO /'EY 
films must take into account a remarkably large diffusion coefficient for eosin in 
solution. A quantitative comparison with the ratio Do/ De of the model is not 
reliable because the model contains a small number of adjustable parameters, 
while the efficiency of other processes such as ZnO precipitation are not taken 
into account. However, it is noticeable that an unexpectedly large diffusion 
coefficient of EY is also obtained experimentally: 3.4 x lO^^cm^s^^ and 1.4 x 
10~^cm^s~^ for O2 and EY, respectively, in the classical deposition condition 
of the films, that is at 70°C and in chloride medium [19]. Consequently, the 
experimental ratio between diffusion coefficients is near 2.4, a much smaller 
value than the ratio between the masses of EY and O2, which is 20. 

The order of magnitude of the growth rate increase in the films of Figs. 8 and 
9, when compared to the films without EA particles, is also consistent with 
experiments: a factor 3 was obtained in Ref. [19], with hydrogen peroxide as 
the oxygen precursor, when Pe/ Po = 10^^, while a factor 5.5 was obtained 
in Ref. [32], with molecular oxygen precursor, when Pe/ Po = 10~^. Moreover, 
the high eosin loading of 2.1% was observed in the latter experiments [32], 
which also agrees with the predictions of our model. 

We also recall that the rules of the model were proposed so that all EA 
branches are connected to the electrode - see the images in Figs. 7, 8 and 
9. At long times, it is also observed that most EA particles belong to some 
of the branches which are exposed above the film surface - see Figs. 8 and 9. 
Accordingly, a small fraction of EA particles is hidden in the deposit. This 
also agrees with the experimental finding that eosin can be almost completely 
removed by a soft chemical treatment of the hybrid films [29,30]. 

Finally, it is interesting to point out that the configurations of EA branches 
are similar to those encountered in a model of diffusion-controlled deposition 
by Burlatsky et al [38], whose growth mechanisms are the same adopted for 
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the E particles in our model. However, while the screening process of the 
larger clusters was the only reason for small clusters to stop growing in their 
model, here this effect is enhanced by the deposition of ZA particles, which 
may eventually cover an EA cluster and suppress its growth (see e. g. Fig. 8c). 
It is also important to notice the similarity between the EA clusters and the 
electrodeposited copper aggregates of Refs. [46,47], which clearly show such 
screening effects. 

From the point of view of competitive growth models, the results of the present 
paper also show nontrivial features. It is usually observed that one of the com- 
peting dynamics is dominant at long times and large lengthscales, such as in 
the widely studied EW-KPZ crossover [48,49,50,51] or in random to correlated 
growth [52] . However, here it was shown that a less trivial association of these 
different dynamics to different chemical species leads to a cooperative behavior 
where features of both dynamics are to some extent preserved. In the present 
model, the competing dynamics are those of MBDLA without bias for the EA 
particles, and of the EW model for the ZA particles, as discussed in Sec. IV. 



5 Conclusion 



We proposed a rate equation model and a microscopic statistical model to rep- 
resent some features of ZnO electrodeposition with organic additives, mainly 
focusing the application to electrodeposition with eosin Y. 

The rate equation model assumes the formation of branches of the additive 
which are extended above the film surface, and is useful to understand the 
basic features of the growth process. Under conditions that the branches grow 
faster than the pure ZnO deposit, it was shown that both structures grow 
with the growth rate of the former. 

With the statistical model, we were able to reproduce several qualitative fea- 
tures of the electrodeposition of ZnO /EY films by assuming that it may be 
viewed as a diffusion-limited aggregation process with an interplay between 
different chemical species in solution and a catalytic effect of the eosin clusters 
for the formation of ZnO. Among the predictions of the model which are in 
qualitative agreement with experimental findings, it is important to mention: 
1) the requirement of particularly large diffusion coefficients of EY in solution, 
so that it is able to improve the growth process with small concentrations; 2) 
the formation of branches of EY connected to the surface and filling the film 
structure without significant large scale inhomogeneities; 3) the increase in 
relative growth rates and the high EY loadings. Consequently, we believe that 
this model incorporated most of the basic features of the ZnO electrodeposi- 
tion with EY, and may be extended in order to provide quantitatively good 
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results or to represent related growth processes. 
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Fig. 1. FESEM cross-sectional view of a ZnO /FiY thin film prepared by electrode- 
position. 
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Fig. 2. Scheme for the rate equation model, with the ZnO film of height h growing 
from the electrode below it, and straight branches of an additive of height H. 




Additive concentration p 



Fig. 3. The solid curve shows the ZnO film growth rate as a function of the additive 
concentration in the rate equation model. The dashed curve is an extension of the 
solid curve for pA < Pc [curve va (pa)]- For higher concentrations of the additive, 
the diffusion of the oxygen precursor becomes limiting, which corresponds to the 
plateau in the growth rate. 
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Fig. 4. (a) Species in the solution of the statistical model are represented by squares 
of sizes equal to 1 lattice unit (O) and 3 lattice units (E). (b) Aggregated species 
(ZA and EA) are represented by squares of the same sizes of the solution. 
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Fig. 5. (a) Rule for aggregation of E particles (E -^ EA), which is possible only upon 
contact with an EA particle, (b) Probabilistic rules for aggregation of O particles 
(O -^ ZA). The aggregation occur with probability 1 in contact with EA. (c) 
Precipitation and diffusion of a ZA particle (indicated by an arrow) immediately 
after its creation. R is the radius to search for the final aggregation point, x indicates 
the point to be chosen in this case, which presents the largest number of neighbors 
ZA. 
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Fig. 6. Section of lateral size 128 of a deposit grown with Dq = 1, De = 1/8, 
po = 0.005, pE = 0.00005, Da = 0.05 and p = 10~^ at i = 3 x 10^ with the 
solution above it. The horizontal line at the bottom represents the electrode. 
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Fig. 7. Initial evolution of a section of lateral size 256 of a deposit grown with 
Do = 1, De = 1, po = 0.02, pE = 0.0002, Da = 0.05 and p = 0.001. Growth times 
are: (a) t = 5 x 10^, (b) t = W^ and (c) t = 1.5 x 10*^. The horizontal line at the 
bottom represents the electrode. 
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Fig. 8. Initial evolution of a section of lateral size 256 of a deposit grown with 
Do = 1, De = 2, po = 0.005, pE = 0.00005, Da = 0.05 and p = 0.25. Growth 
times are: (a) t = lO'^, (b) t = 2 x 10^ and (c) t = 3 x 10*^. The horizontal line at 
the bottom represents the electrode. 
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Fig. 9. Section of lateral size 1024 of a deposit grown with Dp = 1, De = 2, 
po = 0.005, pE = 0.00005, Da = 0.01 and p = 0.25, at t = 2 x 10^. The horizontal 
line at the bottom represents the electrode. 
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